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We have searched for sidereal variations in the rate of antineutrino interactions in the MINOS 
Near Detector. Using antineutrinos produced by the NuMI beam, we find no statistically significant 
sidereal modulation in the rate. When this result is placed in the context of the Standard Model 
Extension theory we are able to place upper limits on the coefficients defining the theory. These 
limits are used in combination with the results from an earlier analysis of MINOS neutrino data to 
further constrain the coefficients. 



PACS numbers: 11.30.Cp,14.60.Pq 



2 



Central to both the Standard Model (SM) and Gen- 
eral Relativity are the principles of Lorentz and CPT in- 
variance. The Standard Model Extension (SME) 0, 0] 
provides a framework for potential Lorentz invariancc 
violation (LV) and CPT invariance violation (CPTV) 
in the SM and suggests such violations could occur at 
the Planck scale, 10 19 GeV. These violations could man- 
ifest themselves at observable energies through several 
unconventional phenomena. One possibility is a poten- 
tial dependence of the neutrino and antineutrino oscil- 
lation probability on the direction of propagation with 
respect to the Sun-centered inertial frame in which the 
SME is formulated [|[. An experiment that has both its 
antineutrino beam and detector fixed on the Earth's sur- 
face could then observe a sidereal variation in the number 
of antineutrinos detected from the beam. 

MINOS is such an experiment jij. It uses Fermilab's 
NuMI neutrino beam [j| and two detectors. The MINOS 
Near Detector (ND) is located 1.04 km from the neutrino 
production target and the Far Detector (FD) is located 
735 km from the production target. The NuMI beam can 
be configured to enhance the muon antineutrino compo- 
nent for high statistics studies using antineutrinos. Both 
detectors are magnetized to approximately 1.4 T, allow- 
ing for the discrimination of /i + produced in charged- 
current (CC) antineutrino interactions from \X~ produced 
in CC neutrino interactions. Because of their different 
baselines, the ND and FD are sensitive to different limits 
of the general SME formulated for the neutrino sector. 
The predicted SME effects for baselines of about 1 km 
are independent of neutrino mass Q , while for long base- 
lines the effects are a perturbation on the standard mass 
oscillation scenario [7j. MINOS has found no statisti- 
cally significant evidence for these effects with neutrinos 
observed in either its ND [§[ or FD @. The high data 
rate in the ND allows us to expand our search to include 
antineutrinos produced by the NuMI beam. 

According to the SME, for short baselines the proba- 
bility that a Vp_ oscillates to flavor D x , where x is e or r, 
over a distance L from its production to its detection due 
to LV and CPTV is given by Q 

Pvn-^v* — L 2 [(C) Sj:i + (AJxp, cos (u®T m ) 

+(A s ) x p, sin (w e T e ) + (B c )xp, cos (2w e T e ) 
+(£ s ) 2p sin(2^T e )] 2 , (1) 

where uj b = 27r/(23' l 56' m 04.0982 s ) is the Earth's sidereal 
frequency, and T® is the local sidereal time of the an- 
tineutrino event. The average value of L is 750 m for 
antineutrinos that are produced by hadron decays in the 
NuMI beam and that interact in the ND. The magnitudes 
of the parameters in Eq. ([TJ depend on the neutrino en- 
ergy, the SME coefficients described below and the di- 
rection of the neutrino propagation in the coordinate 
system fixed on the rotating Earth. The direction vec- 
tors are defined by the colatitude of the NuMI beam line 
X = (90°- latitude) = 42.17973347°, the beam zenith 
angle 6 = 93.2745° defined from the z-axis which points 



up toward the local zenith, and the beam azimuthal angle 
<\> = 203.909° measured counterclockwise from the £-axis 
chosen to lie along the detector's long axis. 

Equation ([1]) for antineutrinos in the ND is identical to 
the oscillation probability equation for neutrinos in the 
ND [8[, with the parameters {A c )xp,, • ■ ■ i (B s )xp replac- 
ing their counterparts (A c ) X fj,, ■ • ■ , {B s ) x ^. The param- 
eter (C) X fi similarly replaces {C) x ^, but does not play a 
role in the sidereal analysis and is not considered further. 

In the SME theory the antineutrino oscillation param- 
eters (A c ) x p,, • ■ • , (B s ) X fi are functions of the coefficients 

( a L)ab ano - ( c i)"f 3- There are 36 of these coefficients: 
the real and imaginary components of (o,l) x , (ai) y , 
(cl) tx , (c l ) ty , (c l ) xx , (c l ) yy , (c l ) xy , (c l ) yz , 
(cl) XZ for f/j, — > v e and Vy, — > v T . Further, these same 36 
coefficients also describe the neutrino oscillation param- 
eters {A c ) xtl , • ■ • , (B s ) xfi . However, the way in which the 
real and imaginary components of the (q.l)" 6 and (cl)^ 
coefficients participate in the (Ac)xp, ■ ■ ■ , (B s ) X fi parame- 
ters is different from the way in which they participate in 
(A c ) xtl , • ■ ■ , {B s ) xlll . The reason for the difference is the 
decomposition of the (a^) and (cl) q/3 coefficients into 
real and imaginary components. For neutrinos 



(a L ) a ab 
and for antineutrinos 



ne{a L T ab + iIm{a L ) a ab 



~Ke{a L )« lb 
ne(c L ) 



\a/3 
ab 



ViTm{a L ) a ah 



ib ■ 



(2) 



(3) 



The subscript "L" in Eq. ([2]) reflects the left-handed na- 
ture of neutrinos while the subscript "i?" in Eq. re- 
flects the right-handed nature of antineutrinos. There is 
a possibility that fortuitous cancellations in the many 
SME coefficients describing neutrino oscillations could 
have masked the sidereal signal for which we were search- 
ing. However, the different dependencies of the param- 
eters for neutrinos and antineutrinos on the SME coef- 
ficients suggest that it is unlikely that a second set of 
fortuitous cancellations would also mask an LV sidereal 
signal for antineutrinos. 

Our primary motivation for this analysis is to explore 
a new window into LV with antineutrinos. Furthermore 
this analysis sheds light on whether cancellations among 
the SME coefficients can affect the results. If MINOS 
is sensitive to sidereal effects resulting from LV in the 
neutrino sector and these effects are being masked by 
accidental cancellations, then this antineutrino analysis 
would find them. On the other hand, if we find no signif- 
icant evidence for a sidereal signal in antineutrinos, we 
can use our results to improve the MINOS upper limits 
on the SME coefficients we previously found with neu- 
trinos since the same coefficients describe both neutrino 
and antineutrino oscillations. 

We applied standard MINOS beam and data quality 
selection [l(| to select beam spills for the analysis. We 
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also applied data quality cuts to remove data where there 
were cooling system problems, magnetic coil problems, or 
an incorrectly configured readout trigger. 

Two independent periods of muon antineutrino data 
taking are combined to comprise the data set for this 
analysis. Table U gives the run dates, number of protons 
incident on the target (POT), and the number of CC 
events remaining in the sample after all selections have 
been made, Ncc- The events were selected following the 
prescription of a previous MINOS analysis [HI]. Stud- 
ies have shown that the mean number of antineutrinos 
per POT in the ND has remained stable to about 1% 
throughout the data taking. Our previous analysis used 
3.54 x 10 6 muon neutrinos observed in the ND Q. 



TABLE I: Antineutrino Data Sample. 



Run Dates 


POT 


Ncc 


Sep09 - MarlO 
NovlO - Janll 


1.67 x 10 2U 
0.98 x 10 20 


637,805 
379,877 


Total 


2.65 x 10™ 


1,017,682 



We used the ratio of the events observed to the num- 
ber of POT recorded as a function of sidereal time as 
the normalized quantity in which to search for sidereal 
variations. We implemented the search for a sidereal sig- 
nal as a blind analysis where we only examined the event 
rate for the data once the analysis procedures were de- 
termined. We used the sidereal time distribution of the 
beam spills and the total number of antineutrino events 
in the data set as inputs to generate 10 4 numerical exper- 
iments that simulated the data set without a sidereal sig- 
nal. We then performed a Fourier analysis on these simu- 
lated experiments to establish the search criteria needed 
to find a sidereal signal. 

We constructed our simulated experiments based on 
the local sidereal time (LST) distribution, T©, of the 
beam spills converted to local sidereal phase (LSP), 
where LSP = mod (TqU)®/2it). To generate this his- 
togram, we converted the time of the extraction magnet 
signal that initiates each spill, as recorded by a GPS unit, 
into LST in standard ways [12J. The GPS time is accu- 
rate to 200 ns [lH and event times were not corrected for 
the time within the 10 /zs spill. We then computed the 
LSP for each beam spill and entered it into a histogram 
with 32 bins ranging from — 1 in LSP. We chose this 
binning because the Fast Fourier Transform (FFT) al- 
gorithm used to look for sidereal variations works most 
efficiently for 2 U bins [14]. Since Eq. (JTJ only puts 
power into the four harmonic terms u}§T§, . . . jAuj^T®, 
we adopted J\f — 5 as the binning that retains these har- 
monic terms while still providing sufficient resolution in 
sidereal time to detect a signal. Each phase bin spans 
0.031 in LSP or 45 minutes in sidereal time. 

To construct the simulated experiments we took each 
spill in the data set one at a time and randomly assigned 
a new LSP for the spill from the LSP distribution of all 
spills. We assigned the number of POT in the spill to 



one histogram in LSP using the newly assigned phase 
and then checked whether any antineutrino events were 
recorded for the spill. If so, we put those events in a 
second histogram using the same LSP. This procedure 
ensured that the correlation between POT and events 
observed in the ND for each spill was retained. By the 
end of the simulation, we have two histograms: one with 
POT as a function of LSP and one with the events as a 
function of LSP. By picking spill times out of the LSP 
distribution for the data, we are assured that both his- 
tograms have their entries distributed properly in LSP. In 
addition, we guaranteed that no sidereal signal is present 
in the simulated experiments since any correlation be- 
tween the data spills is removed. We took the ratio of 
these two histograms to obtain the rate histogram for the 
simulated experiment. 

We next performed an FFT on each simulated rate his- 
togram and computed the power in the four harmonic 
terms (wgT®, . . . , 4w©T0) appearing in the oscillation 
probability, Eq.([T]). Let Si be the power returned by 
the FFT for the first harmonic term sin (uj^Tq) and 
Ci be the power returned for the first harmonic term 
cos (w^T©); similarly define (5*2, C2), . . . ,( £4, C4). Then 
the statistics we used in our search are 

Pi = ^Sl + Cf , ..., P 4 = yfsJ+ C\. (4) 

We added the powers in quadrature to eliminate the effect 
of the arbitrary choice of a zero point in phase at h LST. 
Fig. [T] shows the distribution of px> • • • >P4 f° r the 10 4 
simulated experiments. The distributions for pi, . . . ,p4 




FIG. 1: The distributions for the quadratic sum of powers 
pi , . . . , P4 from the FFT analysis of 10 4 simulated experiments 
without a sidereal signal. The inset shows the distribution 
for pi with a fit to a Rayleigh distribution having a — 0.09 
superposed. 

are quite similar. These distributions are well described 
by a Rayleigh distribution with a = 0.09, showing that 
the powers for the sine and cosine terms of the various 
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harmonics are uncorrelated and normally distributed in 
the experiments. 

Our threshold for signal detection in any harmonic is 
the quadratic power p(FFT) that is greater than 99.7% 
of the entries in its p\ , . . . , P4 histogram. We take these 
signal detection thresholds as the 99.7% confidence level 
(C.L.) for the probability that a measured quadratic sum 
of powers for any harmonic was not drawn from a dis- 
tribution having a sidereal signal. These thresholds are 
0.30, 0.30, 0.29, and 0.31 for pi, . . . ,p<±, respectively and 
we adopt p t h = 0.31 as the overall detection threshold. 

We determined the minimum detectable sidereal mod- 
ulation for this analysis by injecting a sidereal signal of 
the form Asin(w^T), where A is a fraction of the mean 
event rate, into a new set of 10 4 simulated experiments 
and repeating the FFT analysis. We found that every 
experiment gave p\ > 0.31 when A = 0.8% of the mean 
rate. Thus, this analysis is sensitive to sub-percent level 
sidereal variations in the mean event rate. 

Once the threshold for the signal detection was deter- 
mined we performed the FFT analysis using the actual 
data event rate as a function of LSP shown in Fig. [5J The 



3.95 




FIG. 2: The phase diagram of the CC antineutrino event rate 
for the ND data. The mean rate of 3.84 events per 10 15 POT 
is superposed and has x 2 /ndf — 21.7/31. 

results for p±, . . . ,p4 are given in Table [TT1 This table also 
shows the probability, Vp, that the measured power is 
due to a noise fluctuation. is the probability of draw- 
ing a value of p\ , . . . , p± from the parent distribution in 
Fig. []] at least as large as found in the data. As none of 

TABLE II: Results for the pi , . . . , p4 statistics for the data 
shown in Fig. [2] The third column gives the probability, Vt, 
that the measured power is due to a noise fluctuation. 



Statistic 


p(FFT) 


Vt 


Pi 


0.12 


0.42 


P2 


0.17 


0.16 


P3 


0.13 


0.35 


Pi 


0.10 


0.54 



the values p±, . . . , p4 exceed our detection threshold, we 



find no evidence for a sidereal signal in the antineutrino 
data set. 

We investigated the sensitivity of our results to sev- 
eral sources of systematic uncertainties. In the previous 
MINOS analyses HEl, the NuMI target was observed to 
have degraded, causing a drop in the event rate through- 
out the exposure. Because of this degradation, we ex- 
amined how linear changes in the event rate over time 
would affect the determination of the detection thresh- 
olds and found such changes had no effect. The NuMI 
target was replaced between the data taking period of 
the previous analyses and this analysis. The new target 
did not show evidence of degradation during the course 
of its exposure. Given that systematic changes in the 
event rate were shown not to affect the previous results 
and that there is no evidence for such changes in these 
data, this source of systematic uncertainty is negligible. 

Potential differences in the event rate for data taken 
during the solar day compared to the solar night are 
another possible source of systematic uncertainty. We 
looked for these effects by searching for systematic dif- 
ferences in the event rate as a function of solar diurnal 
phase. These rates are consistent with no diurnal varia- 
tions and we conclude that diurnal effects are not mask- 
ing a true sidereal signal in the data. 

There is a known ±1% uncertainty in the recorded 
number of POT per spill [T3| that could introduce a mod- 
ulation that would mask a sidereal signal. We introduced 
random variations of this scale in the number of POT 
recorded from each spill and repeated the FFT analysis. 
We observed no change in the detection threshold due 
to these variations. Moreover, we observed no changes in 
the detection threshold due to long term drifts of the size 
±5% over six months. Thus we conclude that the POT 
counting uncertainties cannot mask a sidereal signal. 

As first pointed out by Compton and Getting fl5j |. 
atmospheric effects can mimic a sidereal modulation if 
there were a solar diurnal modulation in the event rate 
that beats with a yearly modulation. Following the meth- 
ods described in [lq |. we found the amplitude of the po- 
tential faux sidereal modulation would be only 0.5% of 
our minimum detectable modulation and therefore would 
not mask a sidereal signal that MINOS could detect. 

In the absence of a sidereal signal, we can determine 
the 99.7% C.L. upper limits on the SME coefficients 

^ e Oi)a&> Xm ( a i)afc' ne (°L)lb > a ,Sd Xm(c L )"f using the 
MINOS Monte Carlo simulation [10j. In this simulation, 
events are generated by modeling the NuMI beam line, 
including hadron production by the 120 GeV/c protons, 
propagation of the hadrons through the focusing elements 
and 675 m decay pipe to the beam absorber, and the 
calculation of the probability that any neutrinos gener- 
ated traverse the ND. The ND neutrino event simulation 
takes the neutrinos from the NuMI simulation, along with 
weights determined by decay kinematics, and uses this 
information as input into the simulation of the ND. 

We determined the confidence limit for an SME coeffi- 
cient by simulating a set of experiments in which we set 
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all but this one coefficient to zero. For the first simulated 
experiment, we injected a negligible LV signal into the 
simulation and constructed the resulting LSP histogram. 
We calculated the survival probability for each antineu- 
trino from its energy, the distance it travels to the ND 
in the simulation and a value for the magnitude of the 
SME coefficient causing a negligible LV signal. We used 
this simulated LSP histogram to compute p±, . . . ,P4 for 
the experiment. We repeated the simulation 1000 times 
to obtain the average value of each P\,. • . ,Pi statistic 
for the value of the chosen SME coefficient. We then 
increased the value of the SME coefficient and recom- 
puted the average value of each p\ , . . . , p± for a second set 
of experiments. We continued the process of increasing 
the value of the SME coefficient until the largest average 
value of any pi,. . . ,pi crossed the detection threshold of 
0.31. We took this value of the SME coefficient to be its 
99.7% C.L. upper limit. We then computed upper limits 
for the remaining SME coefficients in the same way. 

The 99.7% C.L. upper limit of the SME coefficients 
are given in Table IIIII These limits were cross-checked 
by simulating 1000 experiments for each coefficient in 
the table, where that coefficient was set to the deter- 
mined limit and the rest were set to zero. The distri- 
butions of the pi , . . . , p4 statistics for these experiments 
showed the measured values in Table |n] were excluded 
at more than the 99.7% C.L. This table has the same 



TABLE III: The 99.7% C.L. upper limit on SME Coefficients 
for up — > Vx] {cll)" have units [GeV] and {clT 13 are unit less. 



(a L ) x 3.3 x W™ 

(c L ) TX 1.5 x 10' 

(cl) XX 7.8 x 10" 

(c L ) XY 3.9 x 10" 

(c L ) xz 2.3 x 10" 



(a L ) r 3.3 x 10~ 2U 

(c L ) TY 1.5 x 10 

(CL) 
(CL) 



21 

7.8 x 10~ 21 
2.3 x 10~ 21 



form as the 99.7% C.L. tables in [|, [g. We point out 
that for this analysis, as for the previous ND neutrino 
analysis 0, each limit in this table actually represents 
the 99.7% C.L. upper limit on 4 SME coefficients. For 
(ql) x these are: ^(cil);*,, Im(ai)^, 1Ze(a L )* T , and 



lm(a L )*. Similarly, (a L ) Y , (c L ) 



TX 



• (cl) 



xz 



repre- 
sent limits on 4 SME coefficients (the IZe and 2m parts 
of the coefficients for P M — > v e and P M — > D T ). 

By setting all but one SME coefficient to zero to de- 
termine its confidence limit, our method is based on the 
premise that our null detection does not result from fortu- 
itous cancellations of SME coefficients that hide a signal 
of oscillation terms in Eq. ([1]). Since the number of SME 
coefficients is large, this could be an issue. In fact, we 
raised this issue in 8, 9] when wc determined confidence 
limits based on our null detections with neutrinos. But 
when taken together, the null searches for a sidereal sig- 
nal with both neutrinos and antineutrinos make it clear 
that fortuitous cancellations are quite unlikely. Although 
both neutrino and antineutrino oscillations are described 
by the same SME coefficients, the oscillation parameters 



for neutrinos and antineutrinos have different, nonlinear 
dependencies on them. Both sets of oscillation param- 
eters would independently have to cancel. We conclude 
that our method for determining the limits is sound. 

For the 9 SME coefficients TZe(a L ) a and TZe(c L ) afi for 
the channel — > v T , the limits found in |9j are the 
most sensitive we can determine with our analyses of 
the MINOS neutrino and antineutrino data. For the re- 
maining 27 SME coefficients, however, we can improve 
the limits by combining the results from Q with those in 
TableHn] Let (CL) U be the 99.7% C.L. upper limit on an 
SME coefficient determined in || and {CL) D the 99.7% 
C.L. upper limit determined here. We combine the two 
limits as 

l/(CLf = l/(CL)l + 1/(CL)% 

where (CL) is the combined 99.7% C.L. upper limit. The 
most sensitive upper limits we have determined with the 
MINOS neutrino and antineutrino data are given in Ta- 
blc lIVI As discussed, the way we determine the upper lim- 



TABLE IV: The most sensitive 99.7% C.L. upper limits on 
the SME coefficients determined by MINOS neutrino and an- 
tineutrino data; (a,L)" b have units [GeV] and [ch)" b are unit- 
less. Unless otherwise indicated, the limits were determined 
using ND data. 



Coefficient 


ab limit 


ab limit 


Ke(a L ) x 
Xm{a L )ab 


e/i 2.2 x 10~ 2() 
e/i 2.2 x 10" 20 


fir 1 5.9 x 10 -23 
jut 2.2 x 10~ 20 


1le(a L y ab 
lm(a L ) Y b 


efj, 2.2 x 10" 20 
e/i 2.2 x 10" 20 


fir 1 6.1 x 10" 23 
/ir 2.2 x 10~ 20 


Ke{c L y a x 
l™(c L ) T a x 


e/i 9.0 x 10" 2a 
e/i 9.0 x 10" 23 


fir 1 0.5 x 10" 23 
/tr 9.0 x 10" 23 


lm(c L )l Y 


e/i 9.0 x 10" 2a 
e/i 9.0 x 10" 23 


fir 1 0.5 x 10" 23 
/tr 9.0 x 10" 23 


lm(c L ) XX 


e/i 4.6 x 10" 2i 
e/i 4.6 x 10" 21 


fir } 2.5 x 10" 23 
/tr 4.6 x 10" 21 


Tm(ch) YY 


e/i 4.5 x 10" 21 
e/i 4.5 x 10" 21 


fir 1 2.4 x 10" 23 
/tr 4.5 x 10" 21 


lm(c L ) X b Z 


e/i 1.1 x 10 -21 
e/i 1.1 x 10" 21 


jut' 0.7 x 10" 23 
/tr 1.1 x 10" 21 


ne( CL yj 

Zm(cL) Y b z 


e/i 1.1 x 10 -21 
e/i 1.1 x 10~ 21 


HT 1 0.7 x 10" 23 
/ir 1.1 x 10" 21 


Ke(c L ) xr 
lm(c L ) XY 


e/i 2.2 x 10~ 2i 
e/i 2.2 x 10" 21 


iir 1 1.2 x 10" 23 
/ir 2.2 x 10" 21 



'Determined using FD data Q. 

its does not distinguish between the real and imaginary 
parts of the SME coefficients for the oscillation processes 
— > v e and — > v T . This is reflected in Table ITVl 
We compare the 36 limits in Table lTvl with those deter- 
mined by LSND and IceCube. In [8( , we showed that the 
MINOS upper limits determined with only ND neutrino 
data were already more sensitive than those found by 
LSND [17]. IceCube analyzed their data using the simple 
"vector model" Q for the real components of four SME 
coefficients for — > v T transitions, giving TZe(ai J ) x T , 
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TZe(a L )l T < 1.8 x 1CT 23 GeV and Ke(c L )™ , TZe(c L )^ 
< 3.7 x 1(T 27 18]. The IceCube a L -type limits are a 
factor of 3 lower and the CL-type limits 4 orders of mag- 
nitude lower than the MINOS limits reported here for 
these four coefficients. 

We have presented a search for the Lorentz and CPT 
violating sidereal signal predicted by the SME theory 
with antineutrinos detected in the MINOS Near Detec- 
tor. We found no significant evidence for sidereal varia- 
tions in a blind analysis of the data. Furthermore, the 
effects of systematic uncertainties on these results are not 
significant. When framed in the SME theory these 
results lead to the conclusion that we have detected no 
evidence for Lorentz invariance violation in the antineu- 
trino data set. While the large number of coefficients 
describing the theory could fortuitously cancel a side- 



real signal, the MINOS antineutrino and neutrino results, 
when taken together, suggest that this is improbable. 

We computed upper limits for the 36 SME coefficients 
appropriate to this analysis. We then combined these 
with the upper limits we found in our previous analyses, 
and the results are given in Table IIVI MINOS provides 
the lowest limits for 32 of these coefficients. 
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